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We study a quantum dot Josephson junction inside an Aharonov-Bolim environment. The geom- 
etry is modeled by an Anderson impurity coupled to two directly-linked BCS leads. We illustrate 
that the well-established picture of the low-energy physics being governed by an interplay of two 
distinct (singlet and doublet) phases is still valid for this interferometric setup. The phase bound- 
ary depends, however, non-monotonically on the coupling strength between the superconductors, 
causing the system to exhibit re-entrance behavior and multiple phase transitions. We compute 
the zero-temperature Josephson current and demonstrate that it can become negative in the singlet 
phase by virtue of the Coulomb interaction U . As a starting point, the limit of large superconduct- 
ing energy gaps A = oo is solved analytically. In order to tackle arbitrary A < oo and J7 > 0, we 
employ a truncated functional renormalization group scheme which was previously demonstrated to 
give quantitatively reliable results for the quantum dot Josephson problem. 

PACS numbers: 74.50.-|-r, 75.20.Hr 



I. INTRODUCTION 

The low-energy behavior of quantum dot Josephson 
junctions is governed by an interplay of superconductiv- 
ity and the Kondo effect. The physics emerging from 
the competition of these correlation phenomena was dis- 
cussed decades ago in the context of magnetic impurities 
inside superconducting metals. '^^^ If the Kondo energy 
scale Tk is much larger than the superconducting gap 
A, local magnetic moments are screened by virtue of the 
Kondo effect. This causes Cooper pairs to break, and 
the ground state of the system becomes a Kondo rather 
than a BCS singlet. In the opposite limit Tk ^ A, 
Kondo screening is disturbed due to the superconduct- 
ing gap at the Fermi energy, and the ground state de- 
scribes free magnetic moments. At temperature T = 0, a 
first order level-crossing quantum phase transition from 
a non-magnetic singlet to a degenerate (so-called mag- 
netic) ground state is observed if A/Tr- increases. 

In recent years, the rise of nanotechnology allowed 
for sandwiching quantum dots between superconduct- 
ing electrodes and for measuring the equilibrium Joseph- 
son current as a function of well-controlablc microscopic 
parameters. ^^"'^^ This was the motivation to theoreti- 
cally re-investigate the interplay between superconduc- 
tivity and Kondo physics in the context of (Anderson- 
like) models which feature all parameters (and not only 
A/Tr-) necessary to describe the experimental quantum 
dot Josephson junction. ^^"^^ Both the phase boundary 
between the (Kondo) singlet and the (magnetic) doublet 
phases and the supercurrent were calculated using reli- 
able many-particle methods to tackle the vital Coulomb 
interaction ^.18,22,23 particular, is was demonstrated 
that the critical supercurrent as a function of the quan- 
tum dot energy e can be obtained in good agreement with 
experimental data.^^ 

If a quantum dot is placed in one arm of a non- 



superconducting closed Aharonov-Bohm geometry, sig- 
natures of the Fano effect can be experimentally observed 
in mesoscopic systems. ^''"^^ In addition, the interferomet- 
ric setup allows for extracting physical properties which 
cannot be accessed by measurements on the isolated dot 
(such as the transmission phase). ^^"^^ Both situations 
were investigated theoretically using Anderson-like im- 
purity models as well as appropriate many-body meth- 
ods in order to obtain a physical understanding consistent 
with the observed data.'^'^''^^ In contrast, no interferomet- 
ric experiments on quantum dots within a superconduct- 
ing environment have been performed so far. However, 
in consideration of the rapid progress in nano science it is 
reasonable to assume that experimental (e.g., transport) 
data on such setups will become available fairly soon. 

As mentioned above, both the singlet-doublet phase 
transition and the supercurrent of the quantum dot 
Josephson junction were extensively investigated theo- 
retically, most times employing the single impurity An- 
derson model with BCS source and drain leads. ^^"^'^ In 
contrast, the Aharonov-Bohm situation where both su- 
perconductors are in addition directly linked by a hop- 
ping matrix element td has only been partly investi- 
gated. Zhang used a slave-boson mean-field approach 
to compute the supercurrent for the (singlet) situation 
where the Kondo temperature Tk is larger than the en- 
ergy gap A."^^ However, the author fails to obtain cor- 
rect results in the analytically-solveable non-interacting 
case U ~ Q, rendering his results questionable.'^'^ The 
opposite (doublet) situation with Tk < A was studied 
by Osawa, Kurihara, and Yokoshi.'^'^ They employ, how- 
ever, a Hartree-Fock framework which cannot account 
properly for Kondo correlations, the latter being a vital 
ingredient for the problem at hand. Most surprisingly, 
both works do not at all address the question whether 
the general picture of the existence of singlet and dou- 
blet low-energy states survives if the superconductors are 
connected directly, and, if so, how the ground state ac- 
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as well as the coupling terms 
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FIG. 1: (Color online) The interferometric quantum dot 
Josephson junction considered in this paper. 



tually depends on the physical parameters (particularly 
td) of the system. It is the first aim of this paper to clar- 
ify this issue and to demonstrate that the T = 'phase 
boundary' is affected non-monotonically by a finite cou- 
pling td > 0, causing the system to exhibit re-entrance 
behavior and multiple singlet-doublet transitions. Sec- 
ondly, we present reliable results for the zero-temperature 
Josephson current J for arbitrary system parameters (not 
focussing on a specific regime of A/Tx^^) and particu- 
larly illustrate that J can become negative in the singlet 
phase. Our starting point is the so-called atomic limit 
A = oo which can be treated analytically even in pres- 
ence of finite Coulomb correlations. In order to address 
arbitrary A < oo and U > 0, we employ the functional 
renormalization group (FRG). By comparison with accu- 
rate data obtained from the numerical renormalization 
group framework, it was illustrated that this (after trun- 
cation) approximate method succeeds both qualitatively 
and quantitatively in producing the phase boundary and 
supercurrent for the simple quantum dot Josephson junc- 
tion {td = O).^'' In addition, we will demonstrate that for 
td > the FRG scheme benchmarks excellently against 
the analytic result at A = oo, thereby altogether provid- 
ing a reliable tool to study the problem at hand. 

This paper is organized as follows. In Sec. II, we intro- 
duce the single impurity Anderson model with directly 
linked BCS superconducting leads and compute the as- 
sociated non-interacting impurity Green function. The 
limit A = oo is solved analytically in Sec. Ill A, and 
Sec. IIIB is devoted to a short introduction of the FRG 
framework. The phase boundary of the singlet-doublet 
level-crossing phase transition as well as the Josephson 
current are discussed in Sees. IV and V, respectively. We 
conclude our paper with a short summary (Sec. VI). 



II. MODEL 

In order to describe the geometry depicted in Fig. 1, 
we introduce the standard BCS and Anderson impurity 
Hamiltonian 
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Both Cska and da denote usual fermionic annihilation op- 
erators. The quantum dot is characterized by a gate volt- 
age e and a local Coulomb repulsion U between spin-up 
and spin-down electrons. For simplicity, we assume that 
the left (s = L) and right (s = R) BCS leads have equal 
superconducting energy gaps A while exhibiting a finite 
phase difference cj) = 2(j)L = —20^. Both leads are locally 
coupled to each other and to the quantum dot by (real) 
hopping amplitudes td and i, respectively. 

As a first step to solve the quantum many-particle 
problem associated with Eq. (la, lb), it is useful to com- 
pute the non-interacting Green function of the quan- 
tum dot. This is achieved straight-forwardly using the 
equation-of-motion technique. One obtains 
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Here, Ti denote the Pauli matrices, and As determines 
the lead 'self-energy': 

< tdgs{z)Tzgs{z)n - gs{z)Tz 

1 - ^gs[z)T3gg{z)T3 

with the definition L = R, R = L as well as the implicit 
understanding (here and in the following) that the inverse 
is multiplied from the left. The local Green function gs{z) 
of the isolated BCS leads is given by 
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We have assumed the local density of states p = J^k ^(^^ 
£sk)/N and thus the hybridization energy 
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to be constant, implementing the so-called wide-band 
limit. 

Following the arguments presented in Ref. 23, one can 
show that the Josephson current Jg = i{[H,Ns]) (with 
Ns=L,Fi. being the particle number operator of the left 
and right lead, respectively) can be computed from the 
exact expression (taking h ^ 1 and the electron charge 
e = 1 in the following) 



Js = 2TElniTr [tgsd{i(^) - tdGssiiuj)] 



(6) 
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where T is the temperature of the system, and Qad{z) 
and Qss(.z) denote the interacting dot-lead and lead- lead 
Green function, respectively. The first and second term 
of Eq. (6) can naturally be regarded as the impurity and 
direct contribution to the supercurrent. Employing the 
equation-of-motion technique and generalizing the result- 
ing relations by virtue of the Dyson equation for U 0, 
both Gsd{z) and Gss{z) can be expressed in terms of the 
interacting dot Green function G{z) as follows: 



Gsdiz) ^ Asiz)g{z) , 
Gs-siz) = As{z)g{z)Al{z*) 



s{z)Tzgs{z) (7) 



1 - tlgs{z)T3gs{z)T3 



At r = 0, the right-hand side of Eq. (6) can be evalu- 
ated analytically up to second order in td, leading to the 
sinusoidal law J ~ sin((/)) which describes the ordinary 
Josephson junction.'^* 

Employing a gauge transformation, one can show that 
the supercurrent obtained from Eq. (6) fulfills Jl = 
— Jr = J, provided that G{z) is given exactly.^'^ In the 
present paper, we focus exclusively on the situation of 
fully symmetric superconducting leads (featuring equal 
energy gaps A = = Aji as well as equal hybridization 
strengths F = = 2Tji) and can thus refrain (since 
Jl = —Jr trivially holds) from addressing the issue of 
current conservation within the approximate functional 
RG approach introduced below. 



III. SOLUTION STRATEGIES 

In order to study the quantum many-particle prob- 
lem implicated by the Hamiltonian of Eq. (la, lb), we 
proceed as follows. First, we demonstrate that in the 
limit of large BCS gaps A = cx3 one can derive an an- 
alytic expression for the phase boundary describing the 
singlet-doublet phase transition (Sec. Ill A). In order to 
tackle the Coulomb correlation U at arbitrary A < oo, we 
introduce a truncated functional renormalization group 
scheme (Sec. IIIB). By comparison with numerical RG 
data, it was previously shown that the latter provides an 
accurate tool to calculate the phase boundary as well as 
the supercurrent for the (non-interferometric) quantum 
dot Josephson junction. 



A. Analytic treatment of the limit A = oo 

For td ~ 0, it was previously demonstrated (see 
Rcfs. 22,23) that the limit of large superconducting gaps 
A = cx) allows for an analytic treatment even in presence 
of finite Coulomb correlations U. This is still possible for 
the Aharonov-Bohm situation. Namely, at A = oo the 
non-interacting dot Green function [Eq. (2)] becomes 
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where the effective parameters e and A are given by 
td cos{(j)) + il 
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and id = irptd- Including the interacting part, the prob- 
lem is thus equivalent to diagonalizing the effective two- 
particle Hamiltonian 

This can be achieved straight-forwardly by virtue of a Bo- 
golyubov transformation.^® It turns out that the ground 
state of the system is either non-degenerate (a singlet 
which at sufficiently large U can be thought of featuring 
Kondo screening and broken Cooper pairs) or doubly de- 
generate (a 'magnetic' doublet generally associated with 
a free spin), illustrating that this well-known picture is 
still valid in presence of a finite coupling td > 0- By 
comparison of the corresponding many-particle energies 
one can show that the level crossing and thus the zero- 
temperature 'phase transition' is determined by the im- 
plicit equation 



4e2 



4A^ 
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Since cos{(j)) can become negative for < (/) < tt, the 
right-hand side of Eq. (11) is not necessarily a mono- 
tonic function of the bare parameters e/F and id/F, im- 
mediately indicating re-entrance behavior and multiple 
singlet-doublet phase transitions. This will be discussed 
in detail in Sec. IV. 



B. Functional renormalization group approach 

The functional renormalization group is one implemen- 
tation of Wilson's general RG idea for interacting many- 
particle systems. It starts with introducing an energy 
cutoff A into the non-interacting Green function of the 
system under consideration. Here, we choose a multi- 
plicative infrared cutoff 0(|a;| — A) in Matsubara fre- 
quency space. By taking the derivative of many-particle 
vertex functions (such as the self-energy) with respect to 
the cutoff parameter A, one obtains an infinite hierar- 
chy of flow equations, and subsequent integration from 
A = oo down to the cutoff-free system A = leads 
to an in principle exact solution of the many-particle 
problem. In practice, however, the infinite hierarchy 
needs to be truncated, rendering the FRG an approx- 
imate method. In this paper, we employ a truncation 
scheme that keeps the flow equations for the self-energy 
and the two-particle vertex evaluated at zero external fre- 
quencies. The resulting approximation to both quantities 
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FIG. 2: (Color online) The critical interaction strength f/crit as a (non-monotonic) function of the direct coupling td for 
different BCS gaps A and impurity energies e, altogether characterizing the singlet-doublet level-crossing phase transition 
of the Aharonov-Bohm quantum dot Josephson junction. Solid lines where obtained from the FRG approach introduced in 
Sec. IIIB, dashed lines display the analytic result derived in the limit A = oo [see Eq. (11)]. The phase difference between the 
left and right superconducting leads is given by = 0.2-k, OAn, O.Gtv, and O.Stt (from top to bottom at td = 0). The axis of the 
insets are scaled the same as the axis of the corresponding main part. 



is frequency-independent, contains at least all terms up 
to order U, and can be computed numerically with mi- 
nor effort. It was demonstrated in recent works that this 
truncated FRG scheme successfully describes correlation 
effects (e.g., aspects of Kondo physics) in quantum im- 
purity systems.^^''"' In particular, comparison with nu- 
merical RG reference data showed that both the singlet- 
doublet phase transition and the Josephson current of 
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FIG. 3: (Color online) The same as Fig. 2, but for fixed e = 
and different A/F = 1000, 5, 2, 1, and 0.5 (from bottom 
to top). The phase boundary always resembles the analytic 
form derived in the limit A = oo [Eq. (11)], only the size of 
the doublet phase shrinks monotonously with the BCS energy 
gap. 



a (non-interferometric) quantum dot Josephson junction 
can be computed reliably using this framework. ^'^ 

The FRG flow equations for the diagonal and anoma- 
lous part of the self energy S"^ and associated 
with the Hamiltonian of Eq. (la, lb) can be obtained 
by a slight generalization of the derivation presented in 
Ref. 23. They are given by 



Re 



yA _ 



El 

27r 



(12) 



and the flow equation of the effective interaction reads 



Re 
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We have defined the matrix g'^{iuu) via 



g°{i 



(Si)' 



yA 



(13) 



(14) 



The initial conditions to the coupled differential equa- 
tions (12,13) read = 0, = 0, and U^^°° = 
U, and one can carry out a numerical integration using 
standard Runge-Kutta routines in order to obtain the 
frequency-independent FRG approximation S = T,^^^ 
and Ea = S^^^ to the self-energy. Thereafter, the 
Josephson current can be computed from Eq. (6) and 
the approximate impurity Green function g ~ g^^'^ . 
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FIG. 4: (Color online) (a) FRG results for the critical interaction Uait as an (almost quadratic) function of the impurity energy 
e for A/T — 5, 2, and 0.5 (from bottom to top), (b) Critical phase difference <?!>crit separating the singlet (S) from the doublet 
(D) phase as a function of the direct coupling td. The displayed behavior is similar to analytic result derived at A = oo. 



IV. PHASE BOUNDARY 

A. No direct coupling (td = 0) 

For the simple quantum dot Josephson junction [t^ = 
0) it was previously demonstrated that the boundary 
between the singlet and doublet phases of the system 
is, even though roughly being governed by the ratio 
A/Tk, an explicit function of all parameters of the 
system. ^^'^^'^^^^^ The latter are the Coulomb interaction 
[/, the quantum dot energy e, the superconducting gap 
A, the dot-lead hybridization T, and the phase differ- 
ence (j). According to Eq. (11), the atomic- limit phase 
boundary for = is determined by 



4e^ 



AT^ cos2(0/2). 



(15) 



illustrating that an increase of either U ot (f> drives 
the system towards the doublet phase, whereas a non- 
degenerate ground state is energetically favored the more 
e is shifted away from particle-hole symmetry.^^ Func- 
tional and numerical renormalization group calculations 
showed that the overall size of the doublet regime shrinks 
for A < oo, but all parameter dependencies of the phase 
boundary can still be understood in analogy to the case 
of A = oo.^'^ We will now demonstrate that the same 
holds true for the more complicated case of td ^0. 



B. Aharonov-Bohm situation, A = oo 

In order to understand how a direct link between 
the superconductors affects the boundary of the singlet- 
doublet phase transition at A = oo, it is instructive to 
study the parameter dependence of the critical interac- 
tion strength UcHt for the case of particle-hole symmetry 
e — first. The quantity f/dit can be defined unam- 
biguously, whereas, e.g., a critical coupling strength t"'' 



cannot due to the structure of Eq. (11). For 
6 = TT, one obtains 
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'2\2 





TT 



and 



(16) 



illustrating that there is a fundamental difference be- 
tween both cases. At small 0, the system is driven into 
the doublet phase if td is increased [see, e.g., the curve 
for (f) = 0.27r in the inset to Fig. 2(b)]. In contrast, the 
phase boundary Ucmtitd) depends non-monotonically on 
the direct coupling strength for (f> = it. At small values 
of td, Ucnit increases quadratically, acquires a maximal 
value and finally falls off quadratically for large td [see, 
e.g., the data for (j) — O.Stt in the inset to Fig. 2(b)]. 
Thus, a system which is initially in a doublet state can 
be driven into the singlet phase by increasing the cou- 
pling td at fixed J7, but eventually always re-enters the 
doublet phase. 

Using Eq. (11), one can show that the behavior of the 
phase boundary J7crit(id) for arbitrary <j) is always qual- 
itatively similar to the case of either = or = tt. 
The onset of a non-monotonic dependence on td occurs 
for (j) Si 0.297r, implying that one can expect to ob- 
serve re-entrance behavior even if the phase difference 
cannot be controlled precisely. If the gate voltage is 
tuned away from the point of particle-hole degeneracy, 
the phase boundary exhibits an additional extremum [see 
Figs. 2(a,c)]. An additional minimum occurs for either 
large and arbitrary e 7^ or small (j) and e < 0, whereas 
one observes an additional maximum for small (p and 
e > 0. Since the critical value Ucrit{td = 0) is always 
larger than the asymptote J7crit(^d — > oo), the system 
can exhibit a total of three singlet-doublet phase transi- 
tions if the coupling strength td is varied at fixed U (and 
large (f). 

Due to the fact that the right-hand side of Eq. (11) 
is a horizontally shifted quadratic function of the impu- 
rity energy e [sec FRG data for finite A < 00 in Fig. 4], 
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FIG. 5: (Color online) Josephson current J (in units of Jo = eA/h) as a function of the impurity energy e for constant — 0.5n 
and td/T = 0, 0.15, 0.3, 0.6, and 1.0 (from bottom to top at large e). The results were obtained from the FRG framework. In 
presence of a finite coupling td, J(e) acquires a Fano-like lineshape analogous to the linear-response conductance of the ordinary 
Anderson model, and the non-monotonic phase boundary manifests as repeatedly appearing and disappearing discontinuities. 
In addition, one observes that the Josephson current can become negative in the singlet phase (see, e.g., the lower right panel 
which describes a singlet situation at intermediate and large e). The displayed behavior is generic for arbitrary phase differences 



the above-mentioned re-entrance behavior can also be ob- 
served by changing the impurity energy e while fixing all 
other parameters. Tuning e away from particle-hole sym- 
metry can occasionally drive a system which is initially 
in a singlet state into the doublet and then ultimately 
back into the singlet phase. In contrast, the dependence 
of C/crit on the phase difference (jj is always monotonous. 
At small td, a doublet ground state is favored if (j) is 
increased, whereas the opposite holds for larger td [see 
Fig. 2(a-c)]. For e < and td oo, the system is again 
monotonously driven towards the doublet phase if (j) is 
increased. One can analytically demonstrate that the 
crossover between the regimes of 9^ C/crit ^ is char- 
acterized by values of td where the phase boundary is 
completely independent of (j). 



C. Aharonov-Bohm situation, arbitrary A 

The truncated FRG scheme introduced in Sec. Ill B al- 
lows for computing (an approximation to) the self-energy 
but does not yield the many-particle eigenstates of the 
system under consideration (in contrast, e.g., to the nu- 
merical renormalization group). Within this approach, 
the phase boundary is determined from discontinuities 
in the supercurrent with the understanding (based on 
the analytic treatment of the limit A = cc as well as on 
NRG calculations at td = 0^^) that the groundstate is 
non-degenerate in the limit of small U. 

For the simple quantum dot Josephson junction {td = 



0). comparison with NRG data illustrated that the ap- 
proximate FRG scheme describes the phase boundary as 
well as the supercurrent both qualitatively and quanti- 
tatively at small to intermediate Coulomb correlations 
U <8T, whereas at larger U qualitative features of both 
quantities are still captured correctly.^'^ Comparing FRG 
data for large A with the analytic result of Eq. (11) shows 
that the former is also well-suited to tackle the problem 
at hand. The FRG reproduces all characteristics of the 
phase boundary at > correctly [compare the insets of 
Figs. 2(a,b) with the main parts], only the size of the sin- 
glet phase is slightly overestimated. The latter tendency 
was already observed at = 0. 

FRG calculations at finite A < oo demonstrate that all 
parameter dependencies of the phase boundary are sim- 
ilar to the case of A = oo, only the size of the doublet 
regime shrinks [see Figs. 2, 3 and 4(a) for detailed com- 
parisons of C/crit (id) and C/crit (e), respectively]. This is 
again consistent with results for the simple quantum dot 
Josephson junction {td — 0).^^ Since the FRG scheme, 
however, is approximative in U but the critical inter- 
action strength C/crit becomes large for small A, it is 
reasonable to additionally study the phase boundary in 
terms of a different quantity. It turns out that a critical 
phase difference 0crit can always be defined unambigu- 
ously, and that the behavior of 0crit(id7 U) for arbitrary 
A is similar to the atomic-limit solution [see Fig. 4(b)]. 
One can thus conclude that all parameter dependencies 
of the phase boundary can be understood from Eq. (11), 
only the size of the doublet regime shrinks monotonously 
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FIG. 6: (Color online) Josephson current J (in units of Jo = eA/h) as a function of the phase difference (p for particle-hole 
symmetry e = 0. Note that in the lower right panel the system is in a singlet state for t^j « F and large illustrating that the 
current can become negative in this regime by virtue of the Coulomb interaction. 



for finite A < oo. 



V. JOSEPHSON CURRENT 

In this Section, we present zero-temperature FRG re- 
sults for the equilibrium supercurrcnt J flowing through 
the Josephson junction in presence of a finite phase dif- 
ference (j) between the superconducting leads. According 
to Eq. (6), this current can be interpreted to comprise of 
a 'direct' and an 'impurity' contribution. In contrast to 
the phase boundary, it is not determined solely by the 
dot Green function, rendering it impossible to derive an 
analytic result for J in the limit A = oo. Thus, we fo- 
cus exclusively on discussing FRG data for the Josephson 
current, again recalling that this framework was success- 
fully benchmarked against numerical RG reference data 
for td = 0.^'^ 

In order to discuss how a direct link between the super- 
conducting leads affects the supercurrcnt J, it is instruc- 
tive to recall the simple quantum dot Josephson junction 
(td = 0) first. For small A/Tk, the system is in the sin- 
glet phase for all impurity energies e, and the current J(e) 
exhibits a lineshape which resembles the linear-response 
conductance of the ordinary single impurity Anderson 
model (see the td = - curves of Fig. 5). In the op- 
posite limit, J(e) changes discontinuously at some criti- 
cal value iecrit as the system enters the doublet phase. 



The current becomes negative and almost independent 
of the impurity energy. ^"^ In both cases, the evolution of 
J(e) in presence of a finite link > is Fano-like.'^'^''^^ 
In addition, the non-monotonic dependence of the phase 
boundary on the coupling strength td results in multi- 
ple singlet-doublet phase transitions manifesting as the 
appearance and disappearance of discontinuities of J(e) 
(see Fig. 5). One should particularly note that no mat- 
ter how small U or A, the system will always enter the 
doublet phase in the limit of large hoppings td-, provided 
that the impurity energy is not too large. 

For the simple quantum dot Josephson junction {td — 
0), the supercurrcnt is always positive (negative) in the 
singlet (doublet) phase. Both does no longer necessarily 
hold in presence of a finite coupling td- Whereas it is 
rather intuitive that J can become positive in the dou- 
blet regime due to the additional direct link (having in 
mind the ordinary Josephson junction where two super- 
conductors are coupled by a hopping td and J > holds 
for < < tt), one can most notably also observe a 
negative current in the singlet phase,^^ particularly at 
small BCS energy gaps A (see Fig. 5). It is, however, 
imperative to point out that this is solely caused by the 
Goulomb interaction, and the supercurrcnt at U = 
(where the FRG becomes exact) always remains positive 
in the singlet phase. In contrast, Zhang (Ref. 32) obtains 
a negative singlet current yet in the non-interacting limit, 
rendering these results a priori highly questionable. 
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The Josephson current as a function of the phase differ- 
ence (j> displays the same characteristics as J(e). Multiple 
phase transitions manifest as appearing and disappearing 
discontinuities of J{(p) (see Fig. 6) and can be ultimately 
understood from the functional form of the atomic-limit 
phase boundary Eq. (11). In addition, the current can 
become negative in the singlet phase in presence of both 
a direct coupling td and finite Coulomb correlations. The 
actual form of J{(f>) is rather complicated. It is displayed 
for various parameter sets in Fig. 6. 

VI. CONCLUSIONS AND OUTLOOK 

In this paper, we have investigated a quantum dot 
Josephson junction embedded within an Aharonov-Bohm 
environment. By analytically solving the atomic limit of 
large BCS gaps A = oo, we have shown that the low- 
energy physics of this system is governed by an inter- 
play of two distinct (singlet and doublet) phases in com- 
plete analogy with the non-interferometric case where 
both superconductors are not coupled directly. The 
phase boundary, however, depends non-monotonically 
both on the coupling strength and the quantum dot 



energy e. By carrying out functional rcnormalization 
group calculations (which benchmark excellently against 
the atomic-limit result) at arbitrary A, we have demon- 
strated that the overall size of the doublet regime shrinks 
monotonously with the gap size, but the functional form 
of the phase boundary always remains similar to the an- 
alytic expression derived at A = oo. Thus, even if all 
system parameters cannot be adjusted experimentally in 
a precisely controlled way, one can quite generally expect 
to observe re-entrance behavior within an interferometric 
quantum dot Josephson junction. At finite couplings i^, 
the supercurrent J(e) acquires a Fano-like lineshape anal- 
ogous to the linear-response conductance of the ordinary 
Anderson model. Most importantly, we have shown that 
Coulomb correlations can cause J to become negative in 
the singlet phase. 
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In the discussion of the phase boundary we implicitly as- 
sume that the hybridization strength F is chosen as the 
unit of energy. 

Since our FRG scheme turns out to be quantitatively re- 
liable only up to some intermediate Coulomb interaction 
C//r ~ 8 (where Kondo correlations start to become impor- 
tant), we cannot describe the singlet-doublet phase transi- 
tion in the extreme limit of A/F <C 1 (as this would require 
Tk/T ^ 1), no matter which quantity is used to describe 
the phase boundary. However, for the simple quantum dot 
Josephson junction {td = 0) it was demonstrated that 
generic physics shows up for small values of Tk/A (but 
not neccessarily extremely small Tk itsself).^*'^^ Since we 
do not observe any qualitative changes apart from an over- 
all increase of the size of the singlet phase in going from 
A = cxD to A/F ~ 0.2, it is reasonable to assume that the 
same holds for the problem at hand. 

For the simple quantum dot Josephson junction {td = 0) , it 
was previously observed that in the doublet phase the su- 
percurrent obtained from the FRG is completely indepen- 
dent of both the two-particle interaction strength U and 
the impurity energy e.^^ Even though this is an artifact 
of the approximative approach, numerical RG calculations 
showed that physical properties only vary slightly in this 

22 

regime. 

Some generalized perturbative approach^^ could be em- 
ployed in order to get some intuitive understanding why J 
can become negative in the singlet phase. This is however, 
out of the scope of this paper. 



